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Abstract: Since the early 2000s, many satellite passive microwave brightness temperature (BT) archives, 
such as the Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E) BTs, 
have become the useful resources for assessing the changes in the surface and deep soil moistures over 
both arid and semi-arid regions. In this study, we used a new soil effective temperature (Ty) scheme and the 
archived AMSR-E BTs to estimate surface soil moisture (SM) over the Nagqu region in the central Tibetan 
Plateau, China. The surface and deep soil temperatures required for the calculation of regional-scale Ti 
were obtained from outputs of the Community Land Model version 4.5 (CLM4.5). In situ SM 
measurements at the CEOP-CAMP/Tibet (Coordinated Enhanced Observing Period Asia-Australia 
Monsoon Project on the Tibetan Plateau) experimental sites were used to validate the AMSR-E-based SM 
estimations at regional and single-site scales. Furthermore, the spatial distribution of monthly mean 
surface SM over the Nagqu region was obtained from 16 daytime AMSR-E BT observations in July 2004 
over the Nagqu region. Results revealed that the AMSR-E-based surface SM estimations agreed well with 
the zn situ-based surface SM measurements, with the root mean square error (RMSE) ranging from 0.042 to 
0.066 m3/m3 and the coefficient of determination (R?) ranging from 0.71 to 0.92 during the nighttime and 
daytime. The regional surface soil water state map showed a clear spatial pattern related to the terrain. It 
indicated that the lower surface SM values occurred in the mountainous areas of the northern, mid-western 
and southeastern parts of Nagqu region, while the higher surface SM values appeared in the low elevation 
areas such as the Tongtian River Basin, Namco Lake and bog meadows in the central part of Nagqu region. 
Our analysis also showed that the new T,yscheme does not require special fitting parameters or additional 
assumptions, which simplifies the data requirements for regional-scale applications. This scheme 
combined with the archived satellite passive microwave BT observations can be used to estimate the 
historical surface SM for hydrological process studies over the Tibetan Plateau regions. 
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1 Introduction 


It is well known that soil moisture (SM) is a key variable in hydrological and climatic studies 
(Avissar, 1998; Gutzler et al., 2005; Wagner et al., 2007; Gentine et al., 2011; Assouline, 2013; 
Thomsen et al., 2013). The knowledge about SM distribution makes it easy to understand the 
changes in the land surface hydrological cycle and weather evolution on a large scale, particularly 
in the complex areas such as Tibetan Plateau of China. However, it is very difficult and expensive 
to monitor the large-scale distribution of surface SM on the ground. Many studies have shown 
that satellite passive microwave remote sensing is one of the best available approaches for 
monitoring the land surface SM at a large spatial scale (Wen and Su, 2003; Gao et al., 2006; Kerr 
et al., 2010; Brocca et al., 2011). Its feasibility has been demonstrated extensively at both regional 
and continental scales (Wigneron et al., 2003; Paloscia et al., 2006; De Jeu et al., 2008; Kolassa et 
al., 2016; Xie et al., 2017). 

The estimated spatial-temporal distributions of land surface SM by applying satellite passive 
microwave brightness temperature (BT) data can be effectively evaluated against the in situ 
ground measurements (Njoku and Li, 1999; Koike et al., 2000; Brocca et al., 2011; Kolassa et al., 
2016). So far, an enormous number of BT data at high microwave frequencies, such as the 
Tropical Rainfall Measuring Mission/Microwave Imager (TRMM/TMI) BT observations and the 
Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E) BT data, 
have been archived since the early 2000s. Specifically, the TRMM/TMI BT observations at the 
microwave frequency of 10.65 GHz with a spatial resolution of approximately 60 km and the 
AMSR-E BT observations at microwave frequencies of 6.9, 10.7 and 18.7 GHz can be used for 
the estimation of surface SM. However, the AMSR-E BT data archived by the National Snow and 
Ice Data Center (NSIDC) in the USA have a higher spatial resolution (approximately 0.25°) than 
the TRMM/TMI BT data (Knowles et al., 2006). Therefore, the current archive of AMSR-E BT 
data is more qualified for estimating the surface SM (Paloscia et al., 2006; De Jeu et al., 2008). 

Soil effective temperature (Tep) is one of the basic parameters in the estimation of satellite 
passive microwave remote sensing based surface SM. It is defined as the weighted temperature of 
soil emitting layers and depends on the soil temperature and moisture profiles. The uncertainty of 
Te¢ may ultimately affect the final surface SM product. Sabater et al. (2011) reported that an error 
of 5 K in Te can lead to a 5% error in the surface SM product. The behavior of Tey is difficult to 
understand because of the complex soil temperature and moisture information in the deeper soil 
layers. In previous studies, Tep has been described by considering the entire soil medium as the 
plane stratified dielectric layers (Wilheit, 1978; Schmugge et al., 1986; Njoku and Li, 1999; Lee 
et al., 2002). Because of the heterogeneity of soil layers and the difficulty of obtaining detailed 
soil profile information, this plane stratified Tep theoretical scheme (i.e., Te) is almost impossible 
to use in practical applications. Subsequently, a simplified two-layer scheme with one parameter 
(Choudhury et al., 1982) has been well accepted in microwave remote sensing of surface SM on 
the basis of plane stratified dielectric layers theory scheme (Wilheit, 1978). The scheme was 
further improved by considering wavelength and surface SM (Holmes et al., 2006; Wigneron et 
al., 2008). However, these improved schemes require stringent parameter calibration in the field. 

Recently, Lv et al. (2014) developed a new soil effective temperature scheme (hereinafter 
called Lv’s scheme) for estimating Tep in the soil microwave radiative transfer process at the 
L-band. Lv’s scheme comprehensively considers the impact of wavelength, surface SM, sampling 
depth and soil temperature on Tep. Its key parameters can be calculated directly from soil 
temperature and moisture observations or from land surface model outputs. Lv et al. (2016a) 
evaluated this new scheme together with three existing prevailing soil effective temperature 
schemes, i.e., single-parameter Beta SMAP (Soil Moisture Active Passive) scheme, Choudhury’s 
two-layer scheme (the current SMAP’s scheme), and two-parameter Wigneron scheme (the 
current SMOS’s scheme; SMOS, Soil Moisture and Ocean Salinity). Assessment results showed 
that Lv’s scheme is considered as the best approximation of Wilheit’s integral theory scheme. In 
another recent study, Lv et al. (2016b) verified Lv’s scheme and its corresponding multilayer 
form to provide instruction on how to configure sensor mounting depth and identify its optimal 
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value for in situ SM observations. Based on the Maqu SM network configurations, they found that 
Lv’s scheme can note which site of the Maqu network provides a relatively high accuracy when 
estimating Tz. On average, nearly 20% of the T. signal cannot be captured by the existing Maqu 
network in the current installation configuration. This configuration means 31.1% of the Teg 
signal beneath being unknown and implies significant uncertainty when estimating Tep. The 
observation configuration of 5 and 20 cm soil depths has proven to be the best configuration 
because only 7.4% of the T.# signal beneath is unknown. However, above-mentioned studies (Lv 
et al., 2016a, b) only tested the performance of Lv’s scheme at the L-band and did not showed the 
real advantage of the scheme, i.e., its broad applicability at different microwave frequencies. 
There is, in fact, no band limit in the derivation process of Lv’s scheme. And, because of the 
aforementioned characteristic that its key parameters can be calculated directly from in situ 
observations or land surface model outputs, the easy-to-use advantage of this scheme at different 
microwave frequencies is more obvious. 

A better understanding in the behavior of Tep at different microwave frequencies helps to 
improve the estimation accuracy of surface SM using the existing archived satellite passive 
microwave BT data. Either overestimation or underestimation of Te will result in dryer or wetter 
surface SM estimates. As mentioned above, the existing archived satellite microwave BT 
products are high frequency observational data. Even at the C-band, the accurate estimations of 
Teg still depend on the information of soil moisture and temperature profiles. In the study of Lv et 
al. (2016a), the suitability of Lv’s scheme was mainly explored at the L-band. However, the 
derivation process of Lv’s scheme is not limited at a fixed microwave frequency, and it is 
expected that this scheme can be able to be applied at other microwave frequencies, such as the 
C-band. 

In the past, many studies have been carried out to estimate the surface SM using AMSR-E BT 
data (Njoku and Li, 1999; Koike et al., 2000; Paloscia et al., 2006; De Jeu et al., 2008; Brocca et 
al., 2011; Kolassa et al., 2016). These works mainly concentrated on the low elevation areas and 
provide significant references for identifying the uncertainties in the estimations of surface SM in 
these areas. However, there are few related studies on high elevation areas, especially the Tibetan 
Plateau. Therefore, it is critically important to improve the estimations of surface SM and 
promote the adaptability of estimated methodologies under the complex plateau conditions. 

The aims of the present study are to verify the validity and applicability of Lv’s scheme and to 
estimate the surface SM using Lv’s scheme and AMSR-E BT data in the Nagqu region of the 
Tibetan Plateau, China. It is our hope that this study could provide a reference for accurate 
estimating the surface SM in the Nagqu region and also promote the applicability of Lv’s scheme 
at a large scale in the Tibetan Plateau. 


2 Materials and methods 


2.1 Study area 


The Nagqu region with an average altitude of 4824 m a.s.l. is located in the central part of Tibetan 
Plateau, China (30°-36°N, 90°-95°E; Fig. 1). The region is characterized by a cold and semi-arid 
climate with a mean annual precipitation of 430 mm and an annual mean temperature of —1.0°C 
(1954-2014). The annual frost-free days are 20 d, and the annual sunshine hours are 2500-3000 h. 
The northwestern part is bounded by the Kunlun Mountains. The Tanggula Mountains which run 
roughly northwest-southeast cross the middle part of the region, and the Nyaingéntanglha 
Mountains which run roughly northeast-southwest stretches across the southern part of the region. 
The northwestern, middle and southern parts of the region are mountainous steppes. Permanently 
frozen soil is distributed in the mountainous areas and in the Tuotuo-Tongtian River Basin and 
Nujiang River Basin. Marshlands are widely distributed in the plain areas. Regions to the north 
and south of the Tanggula Mountains are alpine swamps and meadows mingled with cultivated 
lands and shrubs. There are more than 40 lakes in the Nagqu region, with Namco Lake being the 
largest one. 
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Fig. 1 Location of the study area (Nagqu region) in the context of Tibetan Plateau (a) and the distributions of 
the 12 CEOP-CAMP/Tibet SMTMS sites in the study area (b). CEOP-CAMP/Tibet, Coordinated Enhanced 
Observing Period Asia-Australia Monsoon Project on the Tibetan Plateau; SMTMS, Soil Moisture and 
Temperature Measurement System; AMSR-E, Advanced Microwave Scanning Radiometer for the Earth 
Observing System. 


2.2 Data preparation 


In situ surface SM measurements were carried out during 2002-2004 at the CEOP-CAMP/Tibet 
(Coordinated Enhanced Observing Period Asia-Australia Monsoon Project on the Tibetan Plateau) 
SMTMS (Soil Moisture and Temperature Measurement System) sites and in 2008 at the BJ site. 
Surface SM was recorded by the SMTMS along the north-south transect to present its average 
situation in the Nagqu region. The locations of the 12 CEOP-CAMP/Tibet SMTMS sites are 
plotted in Figure 1b. Measurements were taken using a time domain reflectometry (TDR) 
technique at soil depths of 4 and 20 cm every 10 min and then the data were averaged every half 
hour. The in situ measurements closest to the overpass time of AMSR-E BT observations were 
selected to validate the satellite-based surface SM estimates. Detailed description such as name, 
location, land-cover type and observation items of the CEOP-CAMP/Tibet SMTMS sites can be 
found in Ma et al. (2003). It should be noted that the data of four SMTMS sites (D66, TTH, 
WADD and Amdo; Fig. 1b) were not included in this study because of short time series or data 
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missing. Therefore, the area-averaged surface SM data of D105, D110, NODA, MS3478, BJ, 
MS3637, MS3608 and ANNI were used for ground truthing to validate the AMSR-E-based 
surface SM estimations at a regional scale (the 3x9 AMSR-E grid region (red-line zone) shown in 
Fig. 1b). 

In this study, the daily AMSR-E BT dataset archived in the NSIDC with a spatial resolution of 
0.25° was used. The dataset contained two files per day indicating the time of data acquisition, 
which was at 01:30 and 13:30 local time for descending (nighttime) and ascending (daytime) 
overpasses, respectively. The daily AMSR-E BT data at frequencies of 6.9, 10.7 and 18.7 GHz 
covering the study area were downloaded from the NSIDC 
(http://nsidc.org/forms/nsidc-0302_or.html) for the same period as the CEOP-CAMP/Tibet 
experiments (i.e., 2002-2004). Daily AMSR-E BT data for a 21x25 gridded area covering the 
entire Nagqu region were extracted from the AMSR-E global dataset. Daily AMSR-E BT data for 
a 3x9 gridded area covering the selected CEOP-CAMP/Tibet SMTMS sites (i.e., the red-line zone 
shown in Fig. 1b) and for a 1x1 gridded area covering the BJ site were extracted for the validation 
of the AMSR-E-based surface SM estimations at regional (Nagqu region) and single-site (BJ site) 
scales. In addition, the soil texture information from in situ soil samples by laboratory analyses, 
as well as the clay and sand percentages and soil bulk density data from the Harmonized World 
Soil Database (HWSD) v 1.1 were utilized to compute the surface soil dielectric constant (SDC). 
Furthermore, we resampled the HWSD data with 1-km - spatial resolution 
(http://www.iiasa.ac.at/Research/LUC/External-World-soil-database/HTML/index.html?sb=1) to 
0.25° to match the resolution of the AMSR-E BT data. 

Commonly, it is difficult to estimate the T at a regional scale due to a lack of deep soil 
temperature information. In this study, we modeled the weighted mean soil temperature in the 
second soil layer (1.75—4.51 cm; central depth of 2.79 cm) and the third soil layer (4.51—9.06 cm; 
central depth of 6.23 cm) by the Community Land Model version 4.5 (CLM4.5) in the 
Community Earth System Model version 1_2_2 (CESM1_2_2). Then, the weighted mean soil 
temperature was used as the iterative initial value of surface soil temperature (71) corresponding 
to a depth of 4 cm. Similarly, the soil temperature in the fifth soil layer (16.56—28.92 cm, central 
depth of 21.22 cm) simulated by the CLM4.5 was used as the iterative initial value of deep soil 
temperature (72). It should be noted that the atmosphere forcing data for CLM4.5 was CRUNCEP 
(Climatic Research Unit National Center for Environmental Prediction) data at a spatial resolution 
of 0.1°x0.1°. The simulation time-step of CLM4.5 was half an hour. Moreover, the CLM4.5 
simulated soil temperatures which had the nearest time relative to the AMSR-E sensor’s overpass 
time were selected and re-gridded to a spatial resolution of 0.25°x0.25° and then used to 
participate in the surface SM estimations. 


2.3 Methodology 


In this study, the surface SM inversion problem can be solved as searching for the optimal 
solution for the forward model (®) by satisfying the following requirement (Pulliainen et al., 


1993): 


2 

BT, -P ; 

min ap) oad , (1) 
f a 


where, f is the microwave frequency (GHz) and three lowest frequencies in AMSR-E (6.9, 10.7, 
and 18.7 GHz) were used in this study; BT; is the satellite-observed microwave brightness 
temperature at frequency f; Px) is the forward model-simulated BT at frequency f, and opis the 
standard deviation of the measurement noise at frequency f. 

The forward model (®) can be obtained from the so-called w-t model (Ulaby et al., 1986). 
Then, Lv’s scheme (Lv et al., 2014) was used to calculate the Tep. The details are given below. 
2.3.1 Formation and simplification of the forward model (®) 

Considering vegetation as a uniform attenuating medium, the satellite-based BT from a sparsely 
vegetated surface may be expressed as follows (Jackson et al., 2002): 
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where p is the horizontal (H) or vertical (V) polarization; T, and T4 are upward and downward 
radiation temperatures from the atmosphere, respectively; Ta is the atmospheric optical depth, i.e., 
the integral of the atmospheric absorption coefficient profile; Tsx is the sky temperature; es, is the 
effective emissivity of rough soil surface; t is the vegetation optical depth; 7; is the soil 
temperature; T. is the vegetation temperature; and œc is the vegetation single scattering albedo. It 
should be noted that the unit measurement of temperature T is expressed as K. 

Equation 2 is the basic form of the œ-t model, which describes the combined microwave 
radiative transfer process from bare soil and vegetation and contains atmospheric impacts. Drusch 
et al. (2001) found that the values of T, and Tz could be assumed to have little impact on the 
computation of BT, even for higher microwave frequencies such as 19 and 37 GHz. They also 
emphasized that the differences in the atmospheric corrections using different approaches could 
be negligible, so they used the approach proposed by Liebe et al. (1993). In addition, the cosmic 
microwave radiation contribution is very low and thus the Tsx can be ignored. Furthermore, the 
atmospheric contribution is negligible at lower microwave frequencies (f<20 GHz) compared to 
the surface (soil and vegetation) radiation contribution. Therefore, the œ-t model can be 
simplified to the following form (Njoku et al., 2003; Wigneron et al., 2003): 


BT, =el expt.) +T, (1—@,)[l—expCr,)]x[1+ d-e,,)expCz,)]. (3) 


According to the Q/H surface roughness model (Wang and Choudhury, 1981) and the Kirchhoff 
law, the effective emissivity of rough soil surface (esp) can be expressed as follows: 


esp =l- [a g Oi, + Oral exp Ch), (4) 


where p and q are the orthogonal polarizations; h is the roughness height (cm), which is a surface 
roughness parameter and is a function of (ks)? and 0 (where k is the wave number (cm”!), s is the 
surface root mean square height (cm), and 0 is the incidence angle (°)); and Q is the polarization 
mixing factor, which is a constant and is a function of s (cm) and horizontal correlation length (/; 
cm). Variables r, and r,y are the smooth soil surface reflectivities at H and V polarizations, 
respectively and can be computed by the Fresnel reflection equations as follows (Ulaby et al., 
1982): 


sp 


(2) 
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where 0 (°) is the incident angle and e is the surface SDC related to the volumetric surface SM 
(cm?/cm). Dobson mixing dielectric permittivity model can be used to compute the surface SDC 
(Dobson et al., 1985). 

When considering the microwave emission from a vegetated surface at higher frequencies 
(f>20 GHz), a forward physical model @; can be established on the basis of Equations 2, 4, 5 and 
6. Corresponding to lower microwave frequencies (f<20 GHz), a concise forward radiative 
transfer model (@;) can be established with Equations 3—6. Therefore, AMSR-E BT can be 
simulated by this forward model (®;) as follows (Verstraete et al., 1996): 

BT, = D; (X, X25- Xn S1 S2198) + Xio (7) 
where y;is the calculation residual (K; i=1, 2, ..., n) representing the forward model parameters 
related to soil, vegetation and surface roughness; and Si; S2 ..., 5,; denote the complex satellite 
sensor configurations including view angle, polarization and frequency. 

From Equation 1 it can be seen that the surface SM estimation is an iterative calculation 
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process including a set of land surface parameters which can minimize the residual y,. Prior to this, 
an appropriate nonlinear iterative minimization algorithm must be specified for the practical 
estimation process. The determined optimal input parameters are just the expected estimation 
results. 

2.3.2 Parameter optimization of the estimation algorithm 


As mentioned before, a reasonable forward radiative transfer model for bare soil or sparsely 
vegetated surface can be established from Equations 3—6. However, besides the estimated variable, 
there are still four unidentified model parameters related to vegetation and surface roughness. 
These unidentified model parameters are vegetation optical depth (te), vegetation single scattering 
albedo (œc), polarization mixing factor (Q) and roughness height (A). The nonlinearity of complex 
land surface water and heat physical processes and the inherent spatial variability of soil physical 
properties will lead to difficulties in quantifying soil microwave emission processes. Therefore, it 
will be rather difficult to simultaneously calculate and calibrate multiple appropriate parameters 
for a forward model inversion at satellite observation scales. Another inherent difficulty in 
simultaneously calculating and calibrating multiple model parameters is the computational 
instability of the estimation process caused by the introduction of any inappropriate model 
parameter into a forward model. Therefore, efforts should be directed to reduce the model 
parameters to overcome these difficulties. The parameter that can comprehensively reflect the 
effects of vegetation and roughness is the one required. A composite influencing parameter of 
vegetation and roughness is defined as a=exp(—2t.-h). Tian et al. (2012) proposed the 
parameterization form of this equation. By introducing the parameter a, the form of Equation 3 
used in Jackson et al. (2002) can be expressed as follows: 


BT, = {epa +(1-@,)d-a*) [1+ (l-e,, a?) Tug, (8) 


sp 
where Tep is the soil effective temperature (°C) and it can be estimated as follows according to 
Lv’s scheme (Lv et al., 2014): 


Lg =Tll-e* |+Te, (9) 


where Tı and T are the CLM4.5-simulated iterative initial values of surface soil temperature and 
deep soil temperature (°C), respectively; and Bı is the optical thickness (cm) related to the 
wavelength (A) and can be calculated as follows: 


n 


4n E 
A We 
where «’ and e” are the real and imaginary parts of the surface SDC, respectively; and Ax is the 
sampling depth for the first layer of the soil profile (cm). 

Equations 4, 5, 6, 8, 9 and 10 form a self-contained forward land surface microwave emission 
model @;. In this forward model, pending model parameters æ and œc can be calibrated with 
optimization techniques using prior knowledge, such as in situ measurements of surface soil 
texture, temperature and moisture, and no ancillary information is required. The iterative process 
for inverting the nonlinear model Ø; can be implemented by different error/cost function 
minimization algorithms such as the Simplex, Quasi-Newton, and Levenberg-Marquart methods 
(Press et al., 1986). In this study, the Levenberg-Marquart method (Levenberg, 1944; Marquardt, 
1963) was used to estimate surface SM from AMSR-E BT observations at microwave frequencies 
of 6.9, 10.7 and 18.7 GHz. The input data for the estimation processes are satellite-based 
microwave BTs and CLM4.5-simulated surface and deep soil temperatures. 


B, = Axx (10) 


2.4 Data analysis 


In this study, the HWSD soil texture data were processed and resampled to match with the 
resolution of the AMSR-E BT data by ArcGIS 10.4 software 
(http://www.esri.com/software/arcgis/explorer-desktop). The AMSR-E BT data were processed by 
ENVI5.4 software (http://www. harrisgeospatial.com/docs/using_envi_WhatsNew.html). In 
situ-measured soil temperature and moisture data were extracted and processed by the Fortran90 
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programs to match with the AMSR-E BT data. The CLM4.5 simulation data were extracted by the 
NCAR Command Language (NCL; http://www.ncl.ucar.edu/overview.shtml) scripts and then 
processed by the Fortran90 programs. 

In this paper, statistical indices such as root mean square error (RMSE) and coefficient of 
determination (R?) were used to measure the differences between in situ-measured values (Xobs) 
and simulated values given by a estimation method or model (xsim). 


RMSE= [2 (ops: = Ximi , a 1) 
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Generally, the measured sequence and the simulated sequence have different probability 
distributions. The multi-order moments of observations need to be calculated to correct the 
simulated sequence using the probability characteristic of the observed sequence. However, as for 
one random observation variable whose probability distribution is unknown, it is very difficult to 
calculate its multi-order moments. By calculating the Cumulative Distribution Function (CDF) of 
the observation variable, we can easily obtain all statistic characteristics of the observed sequence. 
In probability theory, the CDF of a real-valued random variable (X) is the function given by 
Equation 13: 


R 


F(x)=P(XSx), (13) 
where P(X<x) represents the probability that random variable (X) takes on a value less than or 
equal to x. 

The correction process of the simulated sequence by using the CDF curve of the observed 
sequence can be expressed as follows (Reichle and Koster, 2004): 


CDF ops(X' model, corrected)= CDF sim(Xmodel, uncorrected) (14) 


where CDF sim and CDF ops are the CDF curves of the simulated sequence and observed sequence, 
respectively; and x and x’ denote the uncorrected and corrected values of CLM4.5-simulated soil 
temperature, respectively. The strategy for deviation removal is to match the CDF of the 
simulated sequence to the CDF of the observed sequence. After correction, the deviation between 
the two sequences will be removed mostly. 


3 Results and discussion 


3.1 Correction of CLM4.5-simulated soil temperature 


CLM4.5 can provide simulations of the surface and deep soil temperatures at a regional scale, and 
the simulations will be used as the iterative initial values of surface SM estimations. In this study, 
we used CLM4.5-simulated weighted mean soil temperature in the second soil layer and the third 
soil layer as the surface soil temperature (Tı in Eq. 9). Similarly, CLM4.5-simulated soil 
temperature in the fifth soil layer was used as the deep soil temperature (T2 in Eq. 9). Figures 2a 
and b show the time series of in situ-measured and CLM4.5-simulated soil temperatures in the 
surface layer (4 cm) and deep layer (20 cm) at the BJ site in 2008, respectively. Though the 
seasonal variations of both surface and deep soil temperatures within a whole year were simulated 
well by CLM4.5, there were obvious deviations between the diurnal variation processes of 
CLM4.5-simulated and in situ-measured values. As for the surface soil layer, there were a 
maximum absolute deviation of 15.09 K and RMSE of 3.99 K between CLM4.5-simulated and in 
situ-measured soil temperature values, despite a high linear correlation coefficient of 0.95. 
Similarly, as for the deep soil layer, there were also a maximum absolute deviation of 9.68 K and 
RMSE of 3.38 K between CLM4.5-simulated and in situ-measured soil temperatures, despite a 
high linear correlation coefficient of 0.97. Therefore, this significant bias between 


JOURNAL OF ARID LAND 


CLM4.5-simuated and in situ-measured soil temperature values must be removed. According to 
Section 2.4, CDF-matching method can be used to correct CLM4.5-based simulation time series. 
Figures 2c and d present the actual steps in such a correction corresponding to Figures 2a and b, 
respectively. Here black line denotes the CDF curve of in situ-measured soil temperatures and 
gray line denotes the CDF curve of CLM4.5-simulated soil temperatures at the BJ site. 

Subsequently, we obtained the CDF curves of CLM4.5-simulated surface and deep soil 
temperatures over the Nagqu region in July 2004 (Figs. 2e and f). Meanwhile, we obtained the 
corrected CDF curves of CLM4.5-simulated surface and deep soil temperatures over the Nagqu 
region based on in situ-measured surface and deep soil temperature data collected at all 
CEOP-CAMP/Tibet SMTMS sites for correcting the CLM4.5-based regional distributions of 
surface and deep soil temperatures. Based on CLM4.5 simulations, we obtained the spatial 
distributions of monthly mean surface and deep soil temperatures over the Nagqu region in July 
2004 (Figs. 2g and h). Finally, the CLM4.5-based regional distributions of soil temperatures at 
surface and deep layers were corrected by the corrected CDF curves (Figs. 2i and j). 

According to CLM4.5 simulation, the regional distributions of surface and deep soil 
temperatures over the Nagqu region had the similar spatial patterns. It is clear that the 
northeastern part of the Nagqu region had relatively high surface and deep soil temperatures, 
while the lowest surface and deep soil temperatures occurred in the Nujiang River Basin, Tuotuo 
River Basin and Namco Lake. After the CDF-matching correction, both the CLM4.5-simulated 
surface and deep soil temperatures significantly increased over the Nagqu region (Figs. 2g-j). 
Specifically, the corrected surface and deep soil temperatures at the regional scale had an average 
increase of 5.99 and 6.81 K, respectively, compared to the uncorrected values. Therefore, we 
conclude that CLM4.5 simulations of soil temperatures are usually lower than the actual values 
and should be corrected via in situ measurements. 


3.2 Selection of Lv’s scheme 


The value of Tey is significantly dependent on soil moisture and temperature information from 
deep soil layers. At the L-band, the soil moisture and temperature information from soil layer 
below 10 cm contributed as high as 60% of the Tep for most parts of the globe (Lv et al., 2016a). 
However, the Beta SMAP scheme does not consider this contribution, because it is based on the 
average temperature of the surface (0.0-1.8 cm) and deep (1.8-9.9 cm) soil layers derived from 
the Modern-Era Retrospective Analysis for Research and Applications (MERRA)-Land model. 
Currently, NASA’s SMAP mission has updated its T. calculation method from the Beta SMAP 
scheme to the Choudhury’s two-layer scheme (Choudhury et al., 1982), which is a simplified 
form of the Wilheit’s integral theory scheme (Wilheit, 1978). However, this simplification does 
not give a detailed physical explanation. The scheme parameters are comprised of a series of 
empirical constants and need to be calibrated at different respective bands. Furthermore, soil 
temperatures in the 0-10 and 10-20 cm layers derived from the MERRA-Land model were used 
in this scheme (i.e., Choudhury’s two-layer scheme). Similarly, the parameters of the current 
SMOS Tep scheme (Wigneron’s scheme) depend on specific soil characteristics (Wigneron et al., 
2001) and need to be retrieved empirically through the experiments. And, soil temperatures in the 
0-10 and 29-68 cm layers derived from the MERRA-Land model were used in this scheme. 
Therefore, parameter calibration of the SMAP and SMOS Ty schemes may be required 
depending on the field observation depth. 

Compared with the above-mentioned three 7.4 schemes, the primary advantage of Lv’s scheme 
is that it provides clear physical interpretation of Tep and does not require any estimation of 
empirical parameters. In other words, it does not require parameter calibration in the field. 
Additionally, Lv’s scheme has made it possible to identify the optimal sensor installation 
configuration for the existing field observation networks. It can also provide remote sensing 
researchers with a reasonable relationship between ground-sensor installation depth and 
satellite-observation penetration depth in the validation processes of satellite SM products. The 
soil and vegetation in the Nagqu region are very similar to those in the Maqu region, Tibetan 
Plateau. The observation configuration of 5 cm for surface soil layer and 20 cm for deep soil layer 
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Fig. 2 Time series of in situ-observed and CLM4.5-simulated soil temperatures at the (a) surface (4 cm) and (b) 
deep (20 cm) layers at the BJ site in 2008, and the corresponding CDF (Cumulative Distribution Function) curves 
at the (c) surface and (d) deep layers; CDF curves of CLM4.5-simulated (e) surface and (f) deep soil temperatures 
over the Nagqu region in July 2004 and the corresponding corrected CDF curves based on the measurements at all 
CEOP-CAMP/Tibet SMSTS sites; CLM4.5-simulated spatial distributions of monthly mean surface (g) and deep 
(h) soil temperatures over the Nagqu region and the corresponding corrected regional distributions of (i) surface 
and (j) deep soil temperatures by CDF-matching. CLM4.5, Community Land Model version 4.5. x and x’ denote 
the uncorrected and corrected values, respectively. 
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in the Maqu network has proven to be the best configuration because only 7.4% of the Tep signal 
beneath is unknown (Lv et al., 2016b). Therefore, it is reasonable to select the 4 (surface) and 20 
cm (deep) soil layers in the Nagqu region to calculate Tep using Lv’s scheme. 


3.3 Teg estimations 


The daytime and nighttime series of T. derived from Lv’s scheme (Lv et al., 2014) using in situ 
soil temperature measurements at 4 and 20 cm soil depths at the BJ site in 2008 are shown in 
Figure 3. As for daytime, at the microwave frequency of 6.9 GHz, the estimated Tep time series 
and in situ-measured soil temperature time series at 4 and 20 cm soil depths had the RMSE values 
of 1.85 and 6.23 K, respectively. As for nighttime, these RMSE values were 1.22 and 3.27 K, 
respectively. Generally, the deviation between estimated Tg values and in situ-measured soil 
temperatures increased with the increase of soil depth. The differences between estimated Toy 
values and in situ-measured soil temperatures at nighttime were obviously smaller than those at 
daytime. With the increase of microwave frequency from 6.9 to 10.7 and then to 18.7 GHz, all the 
deviations between estimated T. values and in situ-measured soil temperatures consistently 
increased either for daytime or for nighttime, despite the fact that all the increasing amplitudes of 
these deviations were very small. 

In some previous studies on passive microwave remote sensing (e.g., De Jeu, 2003), due to the 
lack of soil temperature information at the regional scale, high microwave frequency BT data 
were used to linearly fit the regional distribution of Tep, sometimes even considering the infrared 
land surface temperature or vegetation temperature derived from satellites such as AVHRR and 
MODIS as the Tep. Such practices, according to our analysis in this study, will lead to significant 
Teg estimation error. 
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Fig. 3 (a) Daytime and (b) nighttime soil effective temperature (Te) derived from Lv’s scheme (a new soil 
effective temperature scheme developed by Lv et al. (2014)) using the in situ-measured soil temperatures at 4 and 
20 cm depths at the BJ site in 2008 


3.4 Estimation and validation of surface SM 


The surface SM was estimated at both regional (the 3x9 grids area) and single-site (the 1x1 grid 
area covering the BJ site) scales with the proposed algorithm and extracted AMSR-E BT data at 
low microwave frequencies. To validate the satellite-based surface SM, we extracted the in situ 
surface SM measurements which had the closest time relative to the AMSR-E sensor’s overpass 
times from the CEOP-CAMP/Tibet database. The validation results of the surface SM estimations 
at regional and single-site scales are presented in Figure 4. Since the periodic changes in the 
scanning orbit of AMSR-E, the Nagqu region cannot be covered by one complete AMSR-E 
daytime or nighttime BT image for every day. As a result, the dates for nighttime validations were 
not the same as the dates for daytime validations. 

It can be seen from Figure 4 that the RMSE values between AMSR-E-based surface SM 
estimations and in situ-based surface SM measurements in the different validation periods were 
very close to each other. A similar situation was also for R?. During the validation period in 2002, 
the values of RMSE and R? were 0.056 m?/m? and 0.82 for daytime and 0.042 m?/m? and 0.71 for 
nighttime, respectively; whereas during the validation period 2003—2004, the values of RMSE and 
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Fig. 4 Comparisons in time series of in situ-based surface SM (soil moisture) measurements and 
AMSR-E-based surface SM estimations during (a) daytime and (b) nighttime at the single-site scale (over the BJ 
site) in 2008, during (c) daytime from 19 June, 2002 to 31 December, 2002 and (d) nighttime from 21 June, 2002 
to 31 December, 2002 at the regional scale (over the selected CEOP-CAMP/Tibet SMTMS sites), and during (c) 
daytime from 1 September, 2003 to 25 August, 2004 and (f) nighttime from 2 September, 2003 to 26 August, 
2004 at the regional scale. AMSR-E, Advanced Microwave Scanning Radiometer for the Earth Observing System. 


R? were 0.066 m?/m? and 0.92 for daytime and 0.048 m?/m? and 0.90 for nighttime, respectively 
(Figs. 4c-f). Furthermore, during the validation period in 2008, the values of RMSE and R? were 
0.056 m?/m? and 0.78 for daytime and 0.047 m?/m3 and 0.85 for nighttime, respectively. On the 
one hand, as introduced before, there are complicated and changeable geographical environments 
in the Nagqu region. Meanwhile, there are also changeable weather conditions here (Chen et al., 
2012). However, the impact of the complex Asian monsoon climate in the Nagqu region and the 
cloudy weather did not make the atmospheric conditions to affect satellite’s microwave 
monitoring of the ground. In our inversion algorithm, specifically, in Equation 3, the impact of 
the atmosphere on microwave radiation transfer has been ignored. At microwave frequencies of 
6.9, 10.7 and 18.7 GHz, the effect of atmospheric absorption on microwave transfer is often very 
small. Therefore, the simplification about atmospheric absorption in Equation 3 eventually did not 
have any substantial effects on the whole surface SM inversion process. On the other hand, with 
the aid of CLM4.5, we can obtain the detailed temperature distribution information on the deep 
soil layer at the regional scale. Meanwhile, Lv’s scheme can provide us a reliable estimation 
about the whole soil layer microwave emission. These accurate land surface information can 
ensure the stability in surface SM estimations. In fact, the aforementioned validation has proved 
that our retrieval algorithm owns this kind of performance stability in different validation periods. 
Therefore, it can be concluded that this inversion algorithm is basically applicable in the Tibetan 
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Plateau. 

An exact description of the dielectric properties for wet soil is the physical basis of the passive 
microwave estimation of surface SM. Figure 5 shows the time series of surface SDC derived from 
the Dobson soil mixing dielectric model (Dobson et al., 1985) corresponding to three microwave 
frequencies (6.9, 10.7 and 18.7 GHz) using in situ surface SM and soil temperature measurements 
at 4 cm depth at the BJ site in 2008. It can be seen that the dynamic range of surface SDC 
increased with the decrease of microwave frequency. Meanwhile, the real part of surface SDC 
was obviously larger than the imaginary part of SDC. During the whole annual change process, 
the differences in the real part among different microwave frequencies were always significantly 
larger than those in the imaginary part. In addition, the dynamic range of SDC at daytime was 
slightly larger than that at nighttime. By comparing the time series of surface SDC (Figs. 5a and b) 
with those of AMSR-E-estimated and in situ-measured surface SM (Figs. 4a and b), we found that 
the derived surface SDC can adequately reflect the temporal dynamics and seasonal changes of 
surface SM. Therefore, the Dobson soil mixing dielectric model is applicable to the microwave 
SM estimations in the Nagqu region. 
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Fig. 5 (a) Daytime and (b) nighttime series of the real part (rp) and imaginary part (ip) of surface soil dielectric 
constant (SDC) derived from the Dobson soil mixing dielectric model using the in situ surface SM and soil 
temperature measurements at 4 cm depth at the BJ site in 2008 


In the case of all the estimated results in the Nagqu region (see in Fig. 6), the RMSE and R? 
values of AMSR-E-based surface SM estimations and in situ-based surface SM measurements 
were 0.063 m?/m? and 0.90 for daytime and 0.046 m*/m? and 0.83 for nighttime, respectively. 
Therefore, it can be concluded that the AMSR-E-based surface SM estimations agree well with 
the in situ-based surface SM measurements, and the estimated results during nighttime are more 
accurate and stable than those during daytime. 
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Fig. 6 Relationships between AMSR-E-based surface SM estimations and in situ-based surface SM 
measurements for (a) daytime and (b) nighttime during 2002-2004 in the Nagqu region 


However, the comparison between AMSR-E-based surface SM estimations and in situ-based 
surface SM measurements is sometimes questionable due to large differences in spatial 
representativeness (Nicolai-Shaw et al., 2015). The differences in absolute SM between 
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AMSR-E-estimated and in situ-measured values can be related to the numbers of the used 
validation sites, the uniformity of the surface environment, and the spatial size of each satellite 
observation grid. In fact, the basic problem lies in the spatial mismatch between satellite 
observations and ground measurements. In this study, the geographic locations of the field 
investigation sites were carefully chosen so that the AMSR-E-based surface SM estimations could 
be well evaluated against the in situ-based surface SM measurements. At the region scale, a total 
of 27 AMSR-E grids with each area of about 625 km? were involved in the validation process, 
while only 10 validation sites were used in the validation. The situation is worse at the single-site 
scale. Specifically, only one validation site (BJ site) was used to represent the validation area of 
625 km?. Based on above-mentioned analysis, we can say that the deviation between 
AMSR-E-based surface SM estimations and in situ-based surface SM measurements being 0.063 
m/m? for daytime and 0.046 m/m? for nighttime in our validation proves that the 
CEOP-CAMP/Tibet experimental area is one of the most ideal places for satellite SM product 
validation in the world. 

Corresponding to the spatial distributions of monthly mean surface and deep soil temperatures, 
the spatial distribution of monthly mean surface SM over the Nagqu region was obtained from 16 
AMSR-E daytime BT observations on 1, 3, 6, 8, 10, 12, 13, 15, 17, 19, 22, 24, 26, 28, 29 and 31 
July, 2004. The result is shown in Figure 7. 
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Fig.7 Spatial distribution of AMSR-E-estimated monthly mean surface SM over the Nagqu region in July 2004. 
It should be noted that 16 daytime AMSR-E BT observations on 1, 3, 6, 8, 10, 12, 13, 15, 17, 19, 22, 24, 26, 28, 
29 and 31 July, 2004 were used. 


The estimated surface SM values ranged from 0.21 to 0.50 m?/m? in most of the region and the 
values >0.50 m/m? only appeared in the Namco Lake. The Kunlun Mountains in the northern part, 
the Tanggula Mountains in the mid-western part, and the Nyainqgéntanglha Mountains in the 
southeastern part all exhibited low surface SM, with the values ranging from 0.21 to 0.27 m?/m?, 
corresponding to the cold and arid land surfaces. Although the SM estimation algorithm in this 
study did not carefully consider the differences in land-surface types such as glacier/snow cover, 
water body, frozen soil surfaces and the different dielectric behaviors of these land-surface types 
when compared to wet soils, it still successfully and accurately estimated the surface SM in the 
Nagqu region. The Namco Lake and its surrounding area can be clearly identified because of very 
high surface SM. At the spatial scale, as two different segments in the upper reaches of the 
Yangtze River, the Tuotuo River Basin in the high elevation mountainous area exhibited low 
surface SM values, while the Tongtian River Basin in the low elevation area exhibited high 
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surface SM values. Therefore, the AMSR-E BT observations cannot distinguish narrow river 
channels (e.g., the Tuotuo River and the Tongtian River) from their surrounding landscape due to 
the coarse resolution. The lower elevation areas from the mid-eastern part to the southwestern 
part of the Nagqu region mainly exhibited high SM values ranging from 0.27 to 0.50 m?/m?. The 
areas with surface SM values exceeding 0.40 m?/m? may represent the saturated surfaces or water 
bodies in the Nagqu region. Generally, the spatial distribution of monthly mean surface SM 
derived from the AMSR-E BTs clearly reflected the geographical distribution of the land surface 
hydrological status in the Nagqu region. 

The AMSR-E-based surface SM estimations can be further used to understand and analyze the 
historical states or variation trends of surface SM in land surface hydrological processes at large 
spatial scales. The Nagqu region and its adjacent areas are the famous Sanjiangyuan region, 
which is the source region of the Yangtze River, Yellow River, and Lancang River. Therefore, the 
changes in the surface water resources in the region may affect the entire East Asian region. 
Surface SM and regional water resources have a very close hydrological relationship, and thus the 
variations of surface SM may have a direct impact on the regional water resources. Based on the 
surface SM estimations in the Nagqu region, we can obtain the important information on water 
resources in the Tibetan Plateau. 


4 Conclusions 


In this study, a new soil effective temperature scheme (Lv’s scheme) was applied to estimate the 
surface SM over the Nagqu region using the archived AMSR-E BT data at low microwave 
frequencies (f<20 GHz). The in situ-based surface SM measurements were used for the validation 
of the AMSR-E-based surface SM estimations. Generally speaking, Lv’s scheme has a good 
frequency applicability when applied to AMSR-E BT products, despite that it was originally 
designed for the L-band observation mission. The AMSR-E-estimated surface SM agreed well 
with the in situ-measured surface SM both at both regional and single-site scales, demonstrating 
that it is applicable to use Lv’s scheme in the SM estimations over the data-sparse regions such as 
Tibetan Plateau. In the Nagqu region, the derived spatial distribution of monthly mean surface 
SM showed lower SM values in the mountainous areas and higher SM values in the piedmont 
river basins and bog meadow, particularly the areas adjacent to the Namco Lake. The spatial 
distribution of surface SM is closely related to the topographic conditions. Additionally, the 
validation of AMSR-E-based SM estimations showed that the satellite retrieval can capture the 
seasonal variations of surface SM at the regional scale. Although the retrieval method used in our 
study may still need to be improved, the AMSR-E-based SM estimations for different dates 
showed a stable performance, which confirms the potential reduction in attenuation effects of 
atmosphere and clouds on the microwave radiative transfer process. 
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